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We study the phenomenon of macroscopic quantum tunnehng (MQT) in a finite size Josephson 
junction (JJ) with an externaUy apphed magnetic field. As it is well known, the problem of MQT 
in a point-like JJ is reduced to the study of the under-barrier motion of a quantum particle in 
the washboard potential. In the case of a finite size JJ placed in magnetic field, this problem is 
considerably more complex since, besides the phase, the potential itself, depends on space variables. 
We find the general expressions both for the crossover temperature To between thermally activated 
and macroscopic quantum tunneling regimes and the escape time Tcsc . It turns out that in proximity 
of particular values of magnetic field the crossover temperature can vary non-monotonically. 



PACS numbers: 74.50.+r 



I. INTRODUCTION 



The Josephson effecti"^ allows the investigation of fun- 
damental aspects of quantum phenomena such as the 
macroscopic quantum tunneling (MQT)^' which has been 
more recently observed for the first time also in high-Tc 
bi-epitaxial YBCO junctions'*'^. Interesting results have 
been obtained also in various Bi ~ 2212 structures par- 
ticularly referred to "intrinsic" Josephson junctions. 

Usually the phenomenon of macroscopic quantum tun- 
neling is considered in a "point"- like Josephson Junc- 
tions (JJ), i.e. completely neglecting the finiteness of 
the junction size L (some exceptions can be found in 
theoretical and experimental^ papersi»2.). This, zero or- 
der in L, approximation is based on the assumption that 
the junction size is much smaller than all other related 
parameters of the problem, such as the Josephson pen- 

1 /o 

etration depth A,/ = (ftc^/87rejc^cff) and the charac- 
teristic length Lh = £%/£cS, (with £h — [hc/eH] ' 
as the standard quantum mechanical magnetic length)^. 
The effective length l^s depends on the relation between 
the thickness d(i) (z = L, R) of the superconductive elec- 
trodes and the London penetration depth A(i) of the bulk 
superconductor materials, which the electrodes are made 
of. In the limiting cases one can findiS 



-^fl) + ^l'2'l + 



\1) 

d(i) 



(2) 
i(2) 



c^ox, A(i) ^ d(j) 



d{i) < A 



(*) 



(1) 



Thermal fluctuations in Josephson junctions, produce 
a typical " rounding" of the Josephson current branch in 
the TV curvesiiiiSi. Since the pioneering measurements 



of thermal fluctuation phenomena, to obtain such an ef- 
fect, the required condition between thermal energy and 
Josephson coupling energy {Ej '^ T), was usually real- 
ized not by increasing the temperature, rather, by re- 
ducing the value of Ej by applying a proper magnetic 
field. This, even in the case of small J J {L < Xj), has 
significant implications which become of paramount im- 
portance in MQT activation. 

In this context the authors^^ reported the analysis of 
the role of finiteness of the junction's length L obtaining 
the general expression for the crossover temperature Tq 
between thermally activated and MQT regimes for such 
J J. The escape time Tosc was calculated with the expo- 
nential accuracy in the first approximation in [L/ Lh) 
for temperatures T > Tq and for small region below Tq. 
It was demonstrated that the account for the junction's 
size results in the appearance of a strong sensitivity of 
the MQT process on applied magnetic field, making the 
crossover temperature to be non-monotonic function of 
it. Since magnetic field is an easily adjustable parame- 
ter, it can become an important tool in study of such a 
quantum coherent phenomenon without modification of 
other junction parameters. 

In the present article we will proceed and develope the 
study of MQT in a finite size Josephson junction (JJ) 
placed in magnetic field. First we will report the mean 
field solutions for the effective action of such a finite size 
Josephson system and find the values of the effective ac- 
tion at the extremal trajectories. Then, the explicit form 
of phase trajectories close to the extremal ones and corre- 
sponding action functional will be found. This will allow 
us to find the pre-exponential factor in the expression for 



Tcsc of such extended system in a wide temperature re- 
gion, including the crossover point Tg. It turns out that in 
the vicinity of magnetic fields values -ff„ = '^on/ {L£es) 
(<I>o = Trhc/e is the magnetic flux quantum, c is light ve- 
locity), the escape time can vary non-monotonically, a 
new phenomenon which becomes the fingerprint for the 
experimental check of the proposed theory. 

In our analysis we will suppose that L <C Aj but will 
not impose any restrictions on the relation between L 
and Lh- 



II. GENERALITIES 

A. Escape time 

Starting the discussion of the phenomenon of Joseph- 
son current decay in a finite size Josephson junction let 
us recall the substance of this process in a point-like junc- 
tion. 

Let us consider a current biased Josephson tunnel junc- 
tion. In the framework of the capacitively shunted junc- 
tion model it can be represented by the electronic equiv- 
alent circuit, where the resistance of the external circuit 
(-Rcxt), the intrinsic junction capacitance (C) and junc- 
tion resistance (Rt), assumed as a linear ohmic element, 
are connected in parallei^,. The current balance in the 
circuit can be accounted by the following equation 



Ic sin 2(j) - 



where 



V 
Rt 



Vdt 



C 



dV 



(2) 



(3) 



is the relative phase between the two superconductors. 
Equation ([2]) can be rewritten in the form: 



Mc 



d^(t> d(j) dU ((/)) 



dt^ 



where Mc = h^C/e^, -q 



dt 



h{t 



0, 



(4) 



:2i?) 



e 



E J cos 2(1)), 



E.^-^. (5) 



Eq. ([4]) can be considered as the equation of motion with 
friction [r] is its viscosity) of a particle of mass Mq in the 
washboard potential ([5]) (see Fig. [1]). The values of the 
bias current I and the critical current Ic determine the 
slope of the potential [/((/>) and the depth of valleys. 

Let us study the case I < Ic- The minima of potential 
([5]) correspond to the metastable states of the junction. 
We will assume the junction viscosity to be small enough 
not to affect noticeably the particle oscillations in the 
well and its under-barrier motion. At high temperatures 
the thermally activated escape dominates and the result 



U(9) 




FIG. 1: Washboard potential. 



of classical Kramers problem of a particle moving in the 
washboard potential U (x) is valid: 



r(th) 



27r 



exp 



AC/' 



(6) 



where the characteristic plasma frequency LOp is given by 
the interpolation formula 



1 d^U{(l)) 



M, 



c 



/2'KcIc 



1/2 



1- 



1/4 



(7) 



and AC/ = C/max ^ C^min- The extremes of potential C/max 
and C/niin correspond to the local maximum and local 
minimum of the potential energy. 

When the temperature decreases thermal activation 
becomes less and less probable and the process of quan- 
tum tunneling through the barrier becomes important. 
Near some temperature Tq the crossover between Arre- 
nius law and quantum tunneling regime takes place. The 
latter can be considered as the underbarrier motion of 
the particle (instanton propagation). 

At very low temperatures (T <C tUp) the tunneling 
takes place only from the lowest level of the energy spec- 
trum. When /c — / ^ Ic the escape rate can be presented 
in the formM: 



F = 



-(qm) 



1/2 



6ujp /6AC/\ 

VtT \ UJp J 



36AC/ 



dUJr, 



(8) 



One can see that the main difference between Eqs. ^ 
and ([5]) consists only in the temperature dependence of 
the exponent, whereas in the classical case this is the 



activation Boltzmann exponential law with the exponent 
equal to the barrier height divided T, in the case of the 
quantum tunneling temperature in exponent has to be 
substituted by Wp. 

At an arbitrary temperature the combined tunneling 
occurs by the following scheme. First, the particle ex- 
cites in the thermal activation manner and at some mo- 
ment gets up to the energy level En (calculated neglect- 
ing tunneling) and then by means of quantum tunneling 
passes through the barrier. The total tunneling proba- 
bility is determined by the sum over all quantum levels 
of corresponding products of the classical and quantum 
probabilities: 



E 



cxp 



En 



exp{~A{En)}. 



Here 



A{En)= iy/2Mc [Er, 



U{4>)W 



(9) 



(10) 



is the action corresponding to the under-barrier motion 
of the particle with energy En ■ 

One can say that quantum tunneling occurs when the 
particle reaches such a level E^tun at which the proba- 
bility of the direct quantum tunneling through the bar- 
rier becomes larger than the probability of the activation 
jump on the higher energy levels with further tunneling 
through the barrier. Quantitatively this can be formu- 
lated as the condition for the extremum of the exponent 
inEq. ^ 



dA{E) 



dE 



Stun 



1 

T' 



(11) 



Condition pT|) implies that the period of the particle os- 
cillatoin in the inverted potential is equal to 1/T. With 
the growth of temperature the energy _Btun increases and 
when T = w-pjlii it reaches the barrier height. At higher 
temperatures the classical activation escape scheme is re- 
alized. 

The value of the escape time Tcsc of the finite size JJ 
in magnetic field can be determined in the framework 
of the general method of functional integration. In this 
approach the escape rate of the MQT, which is defined 
by the imaginary part of the free energy r'^^ — 2 Im F, 
can be expressed in terms of the partition function of the 
systeroi^: 



F= -TlnZ, 



(12) 



where the latter is defined by the functional integral of 
the S-matrix over all possible trajectories (p (r, r) : 



Z = Dip(r,T)exp{-A 



.^i 



0]) 



(13) 



Here ip (r, r) is the phase difference on the junction at 
the point r and imaginary time r and A [ip (r, r)] is the 
effective action of such extended system. 

The problem of MQT makes sense only within the 
framework of the quasi-classical approximation, i.e. to 
assume that the value of escape time considerably ex- 
ceeds all characteristic time-scales of the internal mo- 
tions. In this case, the imaginary part of the partition 
function is small in comparison with the real one, hence 



2T 



Im / Dip (r, t) exp {~A [ip (r, t)]) 
Re / D(^ (r , t) exp (- A [v:> (r , r)] ) ■ 



(14) 



The imaginary part in the numerator of Eq. (J14p is 
determined by the saddle trajectory of action ip^^^ (r, r) , 
while the real part of the partition function in denom- 
inator has to be calculated at the minimal trajectory 
j^min ^j.^ ^-j _ Pqj. jjjg}j enough temperatures T > Tq, or in 
narrow vicinity of To ( jT — To| <C To), both trajectories 
p^'^^ (r) and 93""" (r) turn out to be time independenti^ii^ 
and the exponential factor in t^J., in analogy with Eq. 
([6]), is determined by the expression 



AA 



CSC 

min 



^ exp [-AALt] , 
= A [p^''' (r)] - A [p" 



(r)] 



(15) 



In the assumption of zero viscosity, one can obtain both 
thermally activation and MQT escaped times (the latter 
even with pre-exponential factor accuracy )i^. For tem- 
peratures above but not too close to Tq it is read as 



■^osc - ^0 



sinh ( 



2TJ 



sin(^ro/r) 



exp{-A<i,} 



(16) 



Hence, in order to find the escape time (Eqs. (fT5l) - p^ ) 
one has to calculate the values of action for trajectories 
<y9^^' and (^™'". The crossover temperature Tq turns to be 
the bifurcation point, below which the time dependent so- 
lution for the saddle point trajectory p^*^^ (r, r) becomes 
more favorable than the static one and the definition of 
escape time presents more sophisticated problems. 



B. Effective action 

The effective action of a Josephson junction placed in 
external magnetic field H as a function of fiowing through 
the junction current /, can be written down in the most 
general case basing on the results o&2J^: 



^['^(r,T)]=- 



1 /-i/sr 

'-' ■/-1/2T 



d, /rf2 I c;^r5^(r,-)\^ ^ 



2e2 V ar 



- 7^'y5(r,T) - ^^ _2 






/•1/2T 

/ dTi{{l-cos[(^(r,T) -(^(r,ri)]}' 

J-1/2T 



trTi2 

"L (t - Ti) ttfl, (r - ri) + cos [yj (r, r) + (^ (r, n)] /3l (r - ti) /3fl (r - n)} + — — r-- 



.1/2T ^ sin2[7rr(r-Ti)] J Sne'^eeS J-1/2T 



he 



\ dr 



[Hx n] 



(17) 



r 



Here C is the junction capacitance, i?N is the tunnel 
resistance, R^h is the shunt resistance, «(!,,/?) (t) and 
l3iL,K) (t) s-re integrated over the energy variable normal 
and anomalous Green functions. The integrals are car- 
ried out over the imaginary time and the junction area. 
The first term corresponds to the kinetic energy of the 
junction. The second and fourth terms describe the con- 
tribution of the potential energy to the action. Let us 
stress that the third term, corresponding to the capac- 
itance renormalization, appears only beyond the quasi- 
classical approximation. The last term in Eq. ([TT]) ac- 
counts for the finite size of the junction and the magnetic 
field contribution to the effective actionji^ 

The natural assumption Tq <C T^ side by side with the 
assumed above condition L <^ Xj allows to considerably 
simplify the general expression Eq. (|17p. what was done 
in the leading approximation in (L/Aj) in Ref.^'^. It 
was shown that in these conditions the third term in Eq. 
(flT)) disappears being reduced to renormalization of the 
capacity C in the first term 



C* ^C 



nh 



duo fdajX fdojA 
2i?N J-00 2n \ div J \ du; J 



(18) 



and changing in the definite way the shape of (p (r, t). 

Variation of Eq. (fT7)) on ip results in getting of the 
quasi-classical equations of motion^", which define the 
extremal trajectories 



SA[cp] 



Sip 



0. 



(19) 



'(r,r) 



Near the extremal trajectory 1^°^^" (r, t) the deviation of 
the function tp (r, r) can be represented in the form of 
expansion by normalized functions ip'^ (r, r) 

^ (r, r) = ^-*^ (r, r) + ^ B^^^'; (r, r) (20) 

which are the eigenfunctions of the equation 



5ip^ 



'^„(r,r)=A>:(r,T). 



(21) 



In this representation the action (jl7[) near the extremal 
trajectory (p'^^^'^ (r, r) is presented by Gaussian type func- 
tional integral over functions (p'l[ (r, r) which can be car- 
ried out analytically. Thus, the problem of definition 



of the value of escape time is reduced to finding of the 
eigenvalues Ajj of the Eq. ((2T|! ^^^. 

We will restrict our consideration within two regions 
of temperatures: a). T > Tq and b) \T - To\ < To. This 
choice is related to the fact that in both these regions one 
can use the functions (p (r, r) in the form Eq. ([20]) with 
time independent ip'^^^'-' (r) . The width of the crossover re- 
gion between thermal activation and macroscopic quan- 
tum tunneling regimes turns to be much smaller than Tq 
and will be estimated below. 



III. EXTREMAL TRAJECTORIES IN STATIC 
APPROXIMATION 



Let us find the explicit form of the extremal trajectory 
j^extr jj-^ static approximation, i.e. in the case when it can 
be considered as time independent. Moreover, in the 
geometry under consideration, when the magnetic field 
is applied along the junction, the phase depends only on 



one coordinate x. Substitution of ip'^^^'^ {x,t) 
in Eq. ( [T9|) leads to the equation 



(^<=^*'' (x) 



.2 av 



2 , -cxtr 



(x) 



1 



+ - sm2ip''''^' (x) 



dx^ ' 2"'" '^ ' ' 2j^S' 

We will look for its solutions in the form 



(p-*'^ {x) 



Vo- -J h (/3 (x) 



(22) 



(23) 



where 939 is a constant, and {ip> (x)) — 0. Corresponding 
boundary conditions are 



dip 
dx 



-L/2 



dip 
dx 



= 0. 



(24) 



= L/2 



One can look for solutions of Eq. (|22p in the 
frameworks of the perturbation theory by the parame- 
ter {L/Xj) . In the first approximation one can rewrite 
Eq. ^ in the form 



^d^^pix) 



I 



2Sjc 



1 
-sin 



^0 



x 



(25) 



This equation with the given above boundary conditions 
is easily solvable: 



Lh . 

—— sm 

Lj 

L 



sm 



L 



2x 

2x 
2a;2 



6i 



H 



LL 



H 



2x L 

Lh Lh 



sin 2iy9o 



cos 2(^0 



(26) 



The phase Lp (x, (/jq) is a periodic function of 930 with 
the period tt. It turns out that some critical value of the 
current Jcr {H) exists such, that for all / < /cr one can 
find for each period two solutions for ip^. The exception 
presents narrow regions of certain magnetic field values, 
which will be found and investigated below. It will be 
shown that for such regions four different solutions for 
</?o exist in the interval ]— 7r/2, 7r/2]. Two solutions cor- 
respond to the minimum of action while the other two 
to its saddle point. When I = I^r pair solutions (one 
corresponding to the minimum, the other to the saddle 
point of the action) confluent in one, while for / > !„ 
the static solution does not exist at all. In the simple 
case of a point-like junction one can find: 



and 



¥'00'" = 2 '"^^^^^^ 



sdl " -^ 



I 

jcS 



I 

jcS 



(27) 



(28) 



When the junction has finite width the analysis is more 
complicated. Nevertheless the smallness L <^ Xj allows 
us to find the relation between ipQ and /, i.e. to find 
j^oxtr ^^^ jy Indeed, integrating Eq. ((22)l over x, one can 
obtain: 

I = jcS / — - sin 2ipo — 



L/2 



L 



2x 
17' 



H 



2^(x) 



(29) 



Substituting in this expression ^p {x) according to Eq. 
()26p and performing integration one finds the equation 
which implicitly relates Lp^ to I : 

2 



Lh 



■ sm 



- — sm 2ifo + 7: 77T— 



^H 



2 V2A 



iff 
L 



sin 4(^0 = —^, 

JcS 



(30) 



where 



L 



L 

SLh . 

sr 

2L 



2L 

l7 



L 

l7 

1 



L 



L 

l7 



(31) 



The critical current Jcr (H) is determined by the value 
/ (^0) , where the point IpQ is solution of the equation 



dli^po) 



dipo 



= 0, 



/cr(-ff)=/('?0). 



(32) 



IV. THE VALUE OF EFFECTIVE ACTION AT 
THE EXTREMAL TRAJECTORIES 

Let us find the value of effective action of the finite size 
J J on the extremal trajectory (p^^^'^ (x). The knowledge 
of A [ip'^'^^'^ (x)] will allow us to find the value of escape 
time with the exponential accuracy in the wide interval 
of temperatures above Tq, or close enough to this bifur- 
cation point (|r - Tol < Tcr^). Further definition of the 
pre-exponential factor will require to perform the func- 
tional integration in Eq. (fT4|) over the trajectories close 
to the extremal one. 

In the case when the phase trajectory p {x, t) = (p (x) 
does not depend on imaginary time (static approxima- 
tion) Eq. (fTT]) is considerably simplified: 



A[^{x)]^ 



eLT 



L/2 r J 

dx <! —-^^^ {x) 

-L/2 



Sjc 




^ • (33) 



^-cos2^(x) + A^^ 



The straightforward integration of this expression with 
the phase 'p{x) determined from the Eqs. ([23]), (| 
results in 



ni hSjc 



8eT 



Lh_ 
2A.7 



2eT 

2 r 



— I sm I — ) cos 2po 



with 



Lh^ 
L 



Lh 
L 



Lh 
L 



L 



L 



cos Aipo 
(34) 



(35) 

This expression already gives in explicit form the value 
of action both for minimal and saddle trajectories (de- 
termined as function of current and magnetic field by 
Eq. jSO])). Let us mention that Eq. dSO]) could be ob- 
tained from Eq. p4p deriving the action and equating 
this derivative with zero: dA (ipo) /dpo = 0. Moreover, 
calcuating the second derivative of the action ((34)) one 
find another useful relation 



d^A (po) _ 2hSjc 
d'^ipo 



eT 

Lh 
2\.j 



^^sin('A)cos2<,o 



L 



h 



cos 4(^0 



(36) 



Note that both Eqs. ^^ and dM]) are valid for any 
value of magnetic field, even in the region where the sec- 
ond (correction) term in Eq. pO|) becomes of the order 
or even larger than the first one. Eq. ([5^ for the value 
of (pQ can be written in the form: 



i5i cos 2(pQ + 2S2 cos Aipo = 0, 



(37) 



where 



Si ~ sin ) and S2 



Solutions of Eq. ((ST]) are 



1 



^H 



2 V 2Aj 



L 



H 



(38) 



The Eqs. ([42])- (gS]) were obtained with the help of Eq. 
([55)1 . Inserting the values of v?^°'" and ipf^ to Eq. ^ 
one can obtain the final expression for the difference of 
actions on the extremal trajectories 



'"'^^' = -^,^^(^2 



(39) 



In the case when |(5i/(52| > 1 the only physically sensi- 
ble solutions are given by the equation 



cos 2^0 



A. 

4^2 



'1 + i 



(40) 



When |(5i/(52| < 1 the solutions corresponding to both 
signs in Eq. ([55]) can be realized. This is the hystere- 
sis domain: the type of solution here depends on the 
prehistory of magnetic field variation. At special points 
Hn = nhcn/ {eL£cg), where L = ttuLh {n = 0,1,2...), 
both states on each period start to be equivalent. 

Let us return to Eq. ([30| . It can be rewritten by means 
of the functions Si {H) and S2 {H) (see Eq. ([55]) as 



AA 



sdl 




which determine in explicit form the exponential factor 
of the escape time (ITSI) . Looking at it one can notice 
the nontrivial oscillatory type dependence of the escape 
time on the value of external magnetic field, which we 
will discuss in the final part of the article. 



Si sin 2(^0 + ^2 sin 4(/3o 



JcS' 



(41) 



This equation can be solved exactly in the algebraic 
functions, but for simplicity we will find its solutions, 
(y5™'" and (/3q'^' in the framework of the perturbation the- 
ory under the assumption \Si\ » 1^21- Simple algebra 
leads to the result: 



and 



1 • M 

— arcsm 

2 \3cSSi 



Si 



I 



jcSSi 



2 XJci'Si 



I 



Si \jcSSi 



, if I/jcSSi > 
ipf\iiI/j,SSi<0 



sdl 



1 ■ M 

— arcsm — 

2 \jcSSi 



^f\iiI/j,SSi>0 
^r\^iI/jcSSi<0 



(42) 



(43) 



V. THE VALUE OF EFFECTIVE ACTION AT 

TRAJECTORIES CLOSE TO THE EXTREMAL 

(PRE-EXPONENTIAL FACTOR) 



The expression ([M)) obtained above allows to deter- 
mine the escape rate with the exponential accuracy, 
what indeed was done above. Determination of the pre- 
exponential factor is more delicate task which requires 
knowledge of the shape of trajectories close to the ex- 
tremal one with further functional integration of the ac- 
tion over them. Now we pass to perform this program. 

In real situation the bifurcation point To lies always 
much below the critical temperature: Tq <C Tc. For tem- 
peratures To < T < Tc, or |To - r| < T^ the general 
expression (|17p can be considerably simplified: 



•^[^<-"i4C/'C'^{^(^^)'-?''""''-^™"'*""" 



AnRshe'^ 



1/2T 



dri 



1/2T 



ip{x,T) -ip{x,T) 



Tl 



f^JcS 2 fdip{x,T) edcsH, 



2 

cxt 



2e 



V dx 



he 



(45) 



Let us find the solutions of Eq. (|21[) in the vicinity of the both time independent extremal trajectories (sur- 



faces). We will look for them in the form: 

(^J; (x, r) = Vr exp (i27rTnT) Xu i^) ■ (46) 

Substitution of this expression in Eqs. ( [T7)) -(PT |) leads 
to equations for x^j [x) : 



+ COS (2^-''- {x))] xl (^) = ^/nXl (x) , (47) 

where we have introduced the Q-factor Q^^ ~ 
2TTTRshC*/h and the parameter C = 2n^T^C*/ (ehjcS) . 
Note, that at the point To one has 



A°±i (To) = 0, 



(48) 



since namely at this point appears the first time- 
dependent solution for the extremal trajectory. 

Let us recall that Eq. (|T7|) is valid for both extremal 
trajectories ip^'^^ and (yS™'". We will look for its solutions 
in the form of perturbation theory series in parameter 
(L/Xj) . There are two sets of corresponding eigenfunc- 
tions. In the zero order approximation they can be odd 
or even in x. 

For even values of fc (fc = 2N, N ^ 0,1, 2...) we have 



xili^) 



/27riV 



- +irn{^) 



L/2 



dxcos 



-L/2 



m 



jZi^)^0- (49) 



For odd values of fc (fc = 2iV + 1, A^ = 0, 1, 2...) the 
eigenfunctions have form 



xirM^)=sm 



7r{2N+l) \ 2Ar+l 



L 



X +ji^:^^ix). 



;^.sm(^i^.),--(.)^0. 

(50) 

The last integrals in Eqs. (|49|) and ([50|) express the 
relations of orthogonality between the first order correc- 
tion 7^„ (x) and corresponding zero approximation solu- 
tion. The functions 7 are supposed to be small enough: 

l7^„(a;)l«l- 
For k ^ one can restrict consideration by the main 

approximation only and get from Eqs. (|48|) . (|49|) . ((50)) the 



following expression for the eigenvalues Aj.„ 



2hj 



— X" — \^ 



(^)' + C("^ + QshH) 



^ff . 


( L \ 


— sm 





L 


K^h) 



i-(¥-fc) J 



cos (2(^o) ■ 
(51) 



For the eigenvalues with fc = (Aj_„) we have to find the 
eigenvalues up to the first correction term in parameter 
{L/\jf. From Eqs. (gS]) , (gH) follows the equation for 

7^„(^): 



g^7^n(:^) 
dx'^ 

Lh . 

-r- sm 



2ipo 



2x 



L 
T^ 



cos 2(^0 = 0. 



(52) 



Its solution, in view of the boundary condition Eq. (I24p . 



is: 



7LW^^ 



Lh 
2A,7 

Lh . 



2x f L 

cos 

Lh \Lh 



sin 2iy9o 



L 
L^ 



2x 
L^ 
2x^ 

^H 



Lh_ 
L 

%L\ 



2x 

L^ 

L 

L^ 

cos2iy9o- (53) 



iFrom Eqs. (HZl) , dH]) , dSS]) one finds the value of A^„ 
with the first correction terms (compare to Eq. ([5T|) ) : 

^A;[„ = C (n' + Qsh\n\) + ^ sin (^j cos 2^0 



+ '.^) -(t^)^°^4^o- 



(54) 



One can notice, that the eigenvalue Ag in the vicinity of 
the saddle point trajectory is negative. This property is 
an obvious consequence of the Eqs. ([5S)l - (|57|) and namely 
this fact results in the appearance of the imaginary part 
of the partition function p^ . 

Substitution of n = 1 to Eq. (|54p gives us the explicit 
definition of the crossover temperature Tiy^: 

C (To) (1 + g,h (To)) + ^ sin (^^ cos 2^30 

^)^(^)cos4^o^0, (55) 

where ipo = ip and it is the function of external current 
/ and magnetic field H (see Eqs. ^- (gH)). For the 
critical value of current, where dl/dfo — 0, Tq = 0. 
This equality follows immediately from equations Eq. 
(|55p. Note that our parameter of perturbation theory 
is {Lh/Xj) and the corrections to eigenfunctions (Eqs. 
(g5)) and ([55| ) are small by this parameter for all values 
of the external current I and magnetic field H. That 
means that the Eq. ([30]) is valid even in the vicinity 
of the points, where sin [L/ Lh) = 0. Important that in 
these regions both terms in the r.h.s. of the Eq. ([50]) 
are of the same order and the nontrivial dependence of 
(fo {I, H) arises, since Eq. ([50]) is equivalent to the fourth 
order equation and in considered region all its coefficients 
turn out to be of the same order. 
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Now one can write down the expression for the effec- 
tive action (|l5l) . valid near both extremal trajectories. 
It is enough to substitute in Eq. (|45l) function ip {x, r) 
in the form of Eq. (PO)) with ip'^ {x, r) defined by Eqs. 
(|46l),(|49l) and ^). Since the eigenvalues A% tend to 
zero as T -^ To for saddle point trajectory we have to 
keep in expansion of the effective action over the coeffi- 
cients {Bj_i} all terms up to the fourth order, keeping 

also the products of the type {{B^Y B'l^, (B°i)^Sf}. 
In result of integration over x and r (see Appendix A) 
the action ([IS]) takes the explicit form as the function of 
coefficients B'l: 



part of the partition function Im Z in Eq. (|14p is deter- 
mined by the integral over trajectories close to the saddle 
point trajectory (/?^'^' (x) (see Eq. (f20| ). Corresponding 
expression for action was already obtained above and it 
is given by Eq. ([56)1 . The functional integral in Eq. ([14]) 
is reduced now to integration over all coefficients i?^. 

Let us start from the integration over the coefficient 
Bq. Since the eigenvalue Aq in vicinity of the saddle point 
trajectory is negative it requires special considerations. 
In order to get the finite answer one has before integra- 
tion to perform analytical continuation Bq — > ib and only 
after that to carry out the integral: 



A[BlBlB\,..Bl.] =A[^,.tr{x)] 



(56) 



Calculation of the coefficients 71^2 is cumbersome but 
straightforward. It is necessary to remember that the 
functions x±n (^) before to be used in Eq. (|45| should 
be normalized. In result one finds 



71 




sin 2iy9o 



(57) 



°° dBg 



/2^ 



exp 



db 

7^ 



lKl{B^Y + 2^,Bt{BlB\) 



exp 



i|A°|62 + 2z7i6(S?S^) 



H 



exp(-^(B?i^%)- 



(59) 



Integration of the coefficients -B^zi; we will yet leave for 
further consideration and now perform that one over the 
coefficients B^j ™^ accordance with the formula 



<PBl , 

27r 



.exp{-AO(sO)' 



-71 



21 



2^ 



[B',B\ 



(60) 



72 



2hj, 



Lh 
2A,7 



T 
2 r 



Lh . L 

sm 

L \Lh 



5k 



cos2ipo- 



— 1 cos 4^0 + 3/x I — 



(58) 



VI. 



OSCILLATIONS OF THE ESCAPE TIME VS 
MAGNETIC FIELD 



Now we are ready to calculate the escape time of the 
"small" J J, which is given by Eq. P^ . The imaginary 



Integrations over the remaining coefficients, besides are 
of the canonical Gaussian type and can be easily per- 
formed. 

What concerns calculation of the real part of the par- 
tition function Re Z it is determined by integration over 
the trajectories passing close to the minimal trajectory 
j^min ^^-j ^-^^ jj-^ ^Yii^ case one can take the action ([56]) 
only with quadratic accuracy over B^ . 

Performing discussed above integrations in real and 
imaginary parts of the partition function (see RefJ-^^)) 
one finds the escape time for high enough temperatures 
T > To, or in the narrow vicinity of Tq (jT — TqI ^ Tq): 



J 




„ ^ exp 

00 27r ^ 



T-i=2Toexp{-A<i,} 

-A;|S?|2-|i30|4(-^2+7?( 



2 



1 



sdl 



^ J mm 



(61) 



with 



^-(2^(5(4 




oc ^ 00 -. 



(62) 



Y.^ 




-^-'-(2A0) ( 



(63) 



Note that the pre-factor m Eq. ((6T|) contains Tq instead 
of T. The eigenvalues {A^}[„iin gdi] are defined by Eqs. 

dnD, HMD- 

The remaining integral over B^ in Eq. (j6ip can be 
expressed in terms of Fresnel integral 



^{x) 



dt exp i-t^ 



(64) 



and one obtains r^g^ in the final form: 



^sdl 




roexp(-AA^°^)yiy2 



A? 



with 



B 



-72 



2\/S 



7? 



exp 



AB 



sdl 



2 

|AS^ 



A" , 

^^2 / sdl 



(65) 



(66) 



Using the explicit Eq. ([M|) for eigenfunctions A'^ one 
can present Eq. (|63p in terms of Euler gamma-function 
r(x): 



1^2 



4;ij,C [r(2-ni(i/))r(2-n2(i/))] 



sdl 



[r(l-ni(H))r(l-n2(H))],„i, 



while the values [^i,2]ga(jdic'["'i,2]min are roots of equation 



C (T) [n2 + Q3h (T) "] + ^ sin T-^ ) cos 2ipo 



Lh 
2\j 



L 



H 



cosiipo — 0, (67) 



written for (po = { (/Jq^'^'^''^ , (/7q"" } accordingly. From Eq. 
(p7)) one can see that the finiteness of L leads to strong 
variation of Y2 and, consequently t~^, as the function of 
magnetic field even for small junction with L -^ Xj. 

All quantities AA^^l^-^,To, t~^ are osciUatory 
functions versus magnetic field (see Eqs.ljTT]), 
(I3ni),(IMl),(Ill,(l53,(IS21))- As the example let us 
consider the behavior of To{H). In first approximation 
by parameter (L/Xj) one can obtain from Eqs. (|¥T|) - ([551) 
following simple expression for the crossover temperature 
To {H) : 



T' 



ejcSRsh 



'S!{H)-[ — 



, 52 [H) 



51 {H) 




with e = QshT = h/ (27ri?shC*) 
Eq. ([38| . Its physical solution 



icS , 

5i and 62 defined by 



n{H) 




AejcSRsh 



Tie 



ISj{H) 



I 

jcS 



/2(g) 



SfiH) 




1/2 



(68) 



evidently oscillates versus magnetic field as it is sketched 
in Fig. [21 Eq. ([68]) can be used until the second term 
in square parenthesis is smaller than the first one. Close 
to the special points {L — ttuLh) this condition can be 
not valid more and in such a case Eas. (l30l) . (|67)) should 
be solved exactly. 

In contrast to Y2, the pre-factor Yi contains only terms 
with Aq (fc y^ 0) and in result depends on magnetic field 
weakly. Nevertheless this dependence turns out to be 
singular due to the logarithmic divergence of the product 
in Eq. (f62|) . Details are presented in Appendix B. 

Eq. ([65]) enable us to estimate the width of the 
crossover region between Arrenius and macroscopic quan- 
tum tunneling regimes. It can be found from condition 



that the argument of the Fresnel function is of the order 
of one: 



{A?}sdi-2x/S. 



VII. FINAL REMARKS 



(69) 



Even a weak external magnetic field can change 
strongly the value of the critical current of the Josephson 
junction of a small size. Note, that some special points 
of external magnetic field appear in the problem under 
consideration. In the vicinity of these points the pertur- 
bation theory fails and the escape time can be calculated 
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only by means of the exact solution of the equations on 
Tq. and ipo- It worth to mention that our general equa- 
tions enable us to consider even such points. 

In the vicinity of the crossover temperature Tq from 
Arrenius law to the quantum tunneling, one can observe 
the strong effect of the finiteness of the junction length 
L even in the pre-exponential factor. 



keep also products of the type {(S?) B'l^, (S" i) Sf }. 
All other terms (containing B'^ with k y^ and n 7^ 0, ±1) 
is enough to take in the second order approximation. In 
this way one writes the value of the effective action ([^5]) 
valid near both extremal trajectories in the form (see 
Rcf.i3) 
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IX. APPENDIX A 

The eigenvalues A'^.^ at the saddle point trajectory 
tend to zero as T — > Tq. In result, when one substitutes 
the function (p (r, r) in the form Eq. (|20p to the action 
Eq. (|45|) and then expands the expression for action in 
Taylor series he has to keep "dangerous" terms, contain- 
ing Bj_i, up to the fourth order. Moreover one has to 





7E/2 71 3k/2 271 |_/|_ 



FIG. 2: Schematic dependence of crossover temperature on 
magnetic field 



A[^(x,t)] = Abcxtr (x)] + IJ2A^\B'J' - ^MllVf f ^Sin(2(pextr (X)) 



B",B"-ixU'-i^B>Sx'o , {B^iYixiy 



IXillllX-il 



k 



llxSlI 



lx?IP 



2hSj, 



T 



^/^ dx 

—- cos [2ipcKti (x, t) 



L/2 



L 



■E 

k 
v2 



^/2 dx 

L/2 



1X^-2 I 



llx?IPIIxMP 



(^°i)'(x°r) 

llxMP 



(x?)'(x^)^ 




(70) 



r 



where ||..|| is the norm of function. Let us carry out the 
integrals entering in the Eq. ([70)l in the explicit form: 



h 



^/2 dx 

— sin (2(/?oxtr [x)) II oiiii 
L IIXillllX-i 



X?X-iXo 



-L/2 



53/2 

Lh 
2\j 



— sm I — ) s\n2ipQ 



L 



sin 4(/5o 



(71) 



and 
h 



'" dx {X\f{x\f 

-./2^'°^^'^-*^^"^^llx?IPIIxMP 



= i{¥^^"(:^'^°^'^" 



1 (Lh 

2 V2A.7 



5k I — — J cos4(y3o + 3^ f — — 



(72) 



with the function /i defined by Eq. p5|) . Using these ex- 
pressions in Eq. ([70|l one can obtain the final expression 
(155)1 for the effective action, valid for trajectories close 
to the extremal ones. Integrals ([7T|) - ([7^ appear in Eq. 
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56)) by means of the functions ((57)1 - (158 



7i 



2nj,s 



Vfh, 



72 



2hjcST 



(73) 



from Eq. (|53)) ). This divergency has the logarithmic 
character and should be cut off at n ~ Tc/Tq and fc ^ 
Aj/^off- As result with the logarithmic accuracy one can 
find: 



X. APPENDIX B 

One can see that expression for the pre-exponential 
factor Yi (|62)) is divergent at {n, k} ^ oo (this follows 

I 






exp ■ 



1 ^ff sin (A ] [cos {2ipl^') - cos (2(pS"")] ^ 



2 L 



iff 



1 



^il-(^7rfc) 

\ ID. T^ / T' 

exp |-^^ Bin (;f ) [cos (2^f ) - cos (2^^)] | E E 

'^ V -ff/ '' fc=l n=-TJTo 



E 



'rrfcAj ^ 



C(^'+Q«n|n|), 



^^^^^ ■ C("' + gs„|n|) 



r 



(74) 



The double sum in the first multiplier in the r.h.s. of this 
expression converges at infinity. 
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